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ABSTRACT 

We develop a formalism to calculate energy-dependent fractional variability (rms) in 
accretion flows. We consider rms spectra resulting from radial dependencies of the 
level of local variability (as expected from propagation of disturbances in accretion 
flows) assuming the constant shape of the spectrum emitted at a given radius. We 
consider the cases when the variability of the flow is either coherent or incoherent 
between different radial zones. As example local emission, we consider blackbody, Wien 
and thermal Comptonization spectra. In addition to numerical results, we present a 
number of analytical formulae for the resulting rms. We also find an analytical formula 
for the disc Wien spectrum, which we find to be a very good approximation to the 
disc blackbody. We compare our results to the rms spectrum observed in an ultrasoft 
state of GRS 1915+105. 
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1 INTRODUCTION 

A very common tool in studying X-ray sources is spectral 
fitting, e.g., with models included in the xspec packet (Ar- 
naud 1996). In this way, a spectrum averaged over some 
period can be decomposed into components due to various 
radiative processes. However, spectral fitting is not unique 
in general, and the same data can often be well fitted with 
more than one spectral model. 

Another way of studying X-ray sources is via their 
power spectra. However, power spectra are usually inte- 
grated over a large range of photon energies, which makes it 
difficult to make inferences regarding the underlying radia- 
tive processes. Also, progress towards physical understand- 
ing of power spectra of, e.g., X-ray binaries, has been rather 
modest. The best present physical description of their power 
spectra appears to be in terms of Lorentzians (e.g., Nowak 
2000; Belloni, Psaltis & van der Klis 2002), Fourier trans- 
forms of damped harmonic oscillators. However, this decom- 
position is rather orthogonal to the decomposition of energy 
spectra into components due to different radiative processes. 

Yet another observable is the energy-dependent vari- 
ance, o~%, or the rms fractional variability, te = cte/Fe, 
where E is the photon energy and Fe is the time-averaged 
flux as a function of E. For a given time bin of the under- 
lying lightcurve and the length of the observation, the vari- 
ance is an integral of the power spectrum over the frequency 
range spanning those two time intervals. Alternatively, one 
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can compute the variance corresponding to a narrower range 
of time scales (frequencies) or to a component in a fit to a 
power spectrum, e.g., the variance corresponding to a QPO. 
Given te and Fe, we can compute the spectrum of the cor- 
responding variable component, 



Fe" = cte = teFe 



G) 



If te is calculated for a range of Fourier frequency, F E al rep- 
resents a frequency-resolved spectrum (Revnivtsev, Gilfanov 
& Churazov 1999) . We note, however, that F™ may not al- 
ways be interpreted as a spectrum from a physical process 
(e.g., for the case of spectral pivoting, see below), while te 
can always be interpreted physically. 

Spectra of rms of astrophysical X-ray sources have of- 
ten been computed (e.g., Mendez et al. 1997; Lin et al. 2000; 
Edelson et al. 2002; Wardzihski et al. 2002; Zdziarski et al. 
2002a, b, 2004, 2005; Vaughan et al. 2003; Gilfanov, Revnivt- 
sev & Molkov 2003; Frontera et al. 2003; Vaughan & Fabian 
2004; Fabian et al. 2004; Gallo et al. 2004; Rodriguez et 
al. 2004). However, interpretation or theoretical models of 
observed rms spectra have so far been relatively rare. E.g., 
Mendez et al. (1997) and Gilfanov et al. (2003) have fitted 
F™ of X-ray binaries by blackbody and Wien spectra, re- 
spectively. Then, Vaughan & Fabian (2004), Fabian et al. 
(2004) and Zdziarski et al. (2005) have modelled observed 
rms spectra due to superposition of different spectral compo- 
nent, each characterized by a different variance (see Section 
Gallo et al. (2004) have shown that the rms spectrum of 
the Seyfert 1H 0707-495 is well reproduced by the ratio of 
the spectra of its high and low flux levels multiplied by a 
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constant. This result can be simply explained as that ratio 
being proportional to the range of fluxes at which the source 
varies at a given energy, which should indeed be proportional 
to its rms. 

As found by Zdziarski et al. (2002b), the dominant 
cause of the variability in the hard state of Cyg X-l, and 
probably some other black hole binaries, on certain time 
scales is the varying flux in the soft seed photons irradiating 
the Comptonizing plasma and much less varying the power 
supplied to the hot plasma. This results in a spectral piv- 
oting with the pivot point in the middle of the spectrum, 
and the rms below the pivot decreasing with energy. Such 
rms spectra decreasing with energy in the X-ray range in the 
hard state have been found, e.g., by Wardzihski et al. (2002) 
and Prontera et al. (2003). Then, Zdziarski et al. (2003) pre- 
sented a formalism to calculate the rms of a pivoting spec- 
trum under various assumptions regarding the distribution 
of the spectral indices and the possible spread of the pivot 
energy. Those authors also show this gives a very good de- 
scription of the long-term rms of Cyg X-f in the hard state. 

An opposite variability pattern has been found in the 
soft state of Cyg X-l. Churazov, Gilfanov & Revnivtsev 
(2001) have found the variability pattern in the RXTE/PCA 
data consists of a constant disc blackbody emission and 
a variable high-energy tail. Zdziarski et al. (2002b) have 
found a very similar pattern in the CGRO/BATSE and 
RXTE/ ASM data, and have interpreted it as resulting from 
a hot plasma with a variable power supply irradiated by an 
approximately constant flux of soft photons. Consistent with 
this picture, they found the rms spectrum of those data in- 
creasing with energy in the X-ray range and saturating at 
an approximately constant level around ~100 keV. 

Another considered cause of energy-dependent rms is 
time evolution of flares. Zycki (2004) have calculated the rms 
spectrum using the flare model of Poutanen & Fabian (1999) 
and showed it results in the rms increasing with energy. 

Here, we further consider the potential of rms spec- 
tra for physical interpretation of X-ray emission of cosmic 
sources and identifying the underlying radiative processes. 
We develop a formalism for calculating rms spectra from 
an accretion flow in which the local level of variability is 
a function of radius, as, e.g., in the propagation models 
of Lyubarskii (1997), Kotov, Churazov & Gilfanov (2001), 
Zycki (2003). In the presented method, we assume that the 
shape of emission from a given radius is constant, and thus 
the changes of the rms with energy are due to superposition 
of emission from different radii. We apply this formalism to 
the cases of blackbody and e-folded power-law local spectra. 
We compare the model rms spectra to one observed from 
the black-hole binary GRS 1915+105. 



2 VARIABILITY FROM DIFFERENT 
SPECTRAL COMPONENTS 

Let us consider emission spectrum averaged over a cer- 
tain period of time consisting of a number of components, 
F — y"V Fi, where Fi is the averege flux in the i-th spectral 
component. The variability of each component is given by 
its standard deviation, <7;. For notational simplicity, we omit 
hereafter the average sign for F and the energy subscript for 
a and F, i.e., F = (Fe), o = cte- 



The total variance, a 2 , has then contributions from all 
of the components, 



(2) 



where the covariance is cry = OiOj in the case of complete 
coherence (same phases), and = for the complete incoher- 
ence. Hereafter, we consider only these two limiting cases, 
in which 



fcoh = / <Ti 



(3) 



respectively. Note that a" coh ^ er inc . 

The relative level of variability in each component is 
n = dijFi. Then, the total rms is 



r C oh = r - j. 



F_ 2 \ - 2 Ff 



(4) 



in the fully coherent and incoherent case, respectively. Note 
that even when the variability of various components is co- 
herent and each n is independent of energy (but different for 
different component), the total r con will change with energy 
due to varying relative contributions of spectral components. 



3 VARIABILITY FROM AN ACCRETION 
FLOW 

We then consider an accretion disc, with the flux per unit 
area, F(R), and the fractional variability, r(R), being func- 
tions of the radius from _Ri n to R ou t ■ The flux from the disc 
per unit angle and energy is 



F = 2tt / dRRF(R). 

- 'Hit, 

If the variability is coherent over that range of R, 



(5) 



a coh = 2n dRRr(R)F(R). (6) 

If the variability is completely incoherent over infinites- 
imally small ranges of the radius, of nc becomes also infinites- 
imally small due to the effect of averaging over a very large 
number of small contributions. On the other hand, the vari- 
ance will be finite when the variability is coherent within 
some finite regions varying approximately independently of 
each other. For example, the variability can be coherent 
within a given coherence length, ARi. Then 



dRRr(R)F(R) 



(7) 



and hereafter the summation covers the range from i?i n to 
.Rout. If ARi is constant (= AR) and small compared to the 
length over which F(R) varies, we have 



oL « 4tt 2 (AJ?) 2 J2 [Rir{Ri)F(Ri)] 2 , 

i=0 



(8) 



where Ri = R ln + (i + 1 /2)AR. Note that although equation 
Q is formally more accurate than the sum of equation JSJ, 
it may not be more accurate in reality given the approxi- 
mate character of the assumption of the sharp boundaries 
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of the independent subregions with constant rms each. We 
also note that the relative rms of equation (JHJ should be 
defined with respect to the flux calculated using the same 
assumptions, i.e., 



F ~ 27tA_R RiF(Ri). 

The sum (S| can in turn be approximated by 



<r inc « An 2 AR / dR[Rr(R)F{R)Y 



(9) 



(10) 



Within this approximation, of nc scales linearly with AR. 
Note that this approximation fails when only few terms con- 
tribute to the sum of equation |(HJ. 

We then consider the case in which the variability is 
coherent over a constant logarithmic range of radius, Aln_R. 
This is likely the case for a self-similar flow, which properties 
are determined by the local size scale. Then equations 0, 
© and 

become 



47t" 



2 E 



din RR 2 r(R)F(R) 



Aln Ri 



47t 2 (Aln R) 2 [R.r(R t )F(R t )] '' 



An 2 Aln R 



i=0 



In R ia 



din R[R 2 r(R)F(R)] 2 , 



(11) 
(12) 

(13) 



where Ri = R[ n exp[(i + 1 /2)Aln R]. In the approximation of 
equation 1131 . a 2 scales linearly with Ami?, which approxi- 
mation fails, however, when only a few regions dominate the 
variability. Equation © is now replaced by 



F « 2n Aln R^2r 2 F(R,) 



(14) 



For comparison, we also calculate the coherent variability 
using the same approximations as for equations 11211 and 
1141 . which yields 



a coh « 27tAln R^2R 2 r(Ri)F(Ri 



(15) 



4 RADIAL DEPENDENCIES OF SPECTRA 
AND THE VARIABILITY AMPLITUDE 

Now we consider some examples of the radiated spectra and 
the radial dependencies of both the spectra and the local 
rms, r(R). For an emitted spectrum, we assume its local 
power is equal to the local dissipation rate, i.e., neglect ad- 
vection. We consider two dimensionless forms, D(R), of the 
dissipation rate, 



D{R) = 
D(R) = 



(R/R in )- :i , 

7 7 [1 - (iW-R) 1/2 ] 

6«(je/ft»)», 



(16) 
(17) 



where Ri n is the inner disc radius. The first dependence cor- 
responds to a Newtonian disc torqued at the inner bound- 
ary (as assumed, in particular, in the commonly-used disc 
blackbody model, diskbb in xspec, Mitsuda et al. 1984). It 
corresponds, e.g., to either R[ n ^> the minimum stable orbit 



(which is between QGM/c 2 and ~ lGM/c 2 from nonrotat- 
ing to the maximally rotating black hole) or the magnetic 
stresses important near the minumum stable orbit (e.g., 
Merloni & Fabian 2003 and references therein). The second 
dependence corresponds to a nonrotating black hole with the 
standard zero-torque inner boundary and Ri n — 6GM/c 2 
(Shakura & Sunyaev 1973). Both expressions are normal- 
ized to the maximum of D = 1, reached at R/R ln — 1 and 
(7/6) 2 , respectively. 

We then consider two simple local spectral forms, a 
blackbody and a power law with both low and high-energy 
cutoffs, which represents a rough approximation to spectra 
from thermal Comptonization. They correspond, e.g., to the 
radiative processes dominant in the soft and hard spectral 
states of black hole binaries, respectively (e.g., Zdziarski & 
Gierliriski 2004). 

The blackbody spectrum is 



F(R) 



AE' 



(18) 



exp[E/kT{R)] - 1' 

where A is the normalization constant, k is the Boltzmann 
constant, and 



T{R) = T max \D{R)] 



1/4 



(19) 



is the local temperature expressed via its maximum, T max . 
We also consider the Wien limit of the blackbody spectrum, 
see Appendix A. 

Then we consider a power law with an exponential cut- 
off, which roughly approximates emission due to unsatu- 
rated Comptonization from a hot plasma. As in the black- 
body case, we normalize the spectrum to the local dissipa- 
tion rate, 



F{R) 



AD(R)E- a( - R) exp(E/E c ) 
V[l-a(R),E Q /E c 



E, 



l-a(R) 



E> E, 



o. 



(20) 



and F = for E < Eq, where F is the incomplete 
Gamma function \T(a,x) — > T(a), x <JC 1, T(a,x) — > 
exp(— x), x ^> 1] a(R) is the local photon index, E c is 
the e- folding energy, the dimensional constant, A, depends 
on the central mass and the accretion rate, and Eq is the 
minimum energy in the spectrum, roughly corresponding to 
the average energy of the seed photons undergoing Comp- 
tonization. 

We then consider a model with a hot inner flow sur- 
rounded by a cold disc, providing the seed photons (e.g., 
Zdziarski et al. 2002). Then the flux in the seed photons de- 
creases with the decreasing radius and the locally-emitted 
spectrum hardens. This is confirmed by Fourier-resolved 
spectroscopy, showing, e.g., the spectrum of Cyg X-l in the 
hard state hardens with the increasing Fourier frequency 
(Revnivtsev et al. 1999). Since the Fourier-resolved spec- 
trum has the slope depending on the logarithm of the Fourier 
frequency, we assume here a(R) to be linear in the logarithm 
of radius, 



a(R) = ct-in + 



ln(R/R in ) 

111 



- am), 



(21) 



where the spectral index changes from a out at the outer 
radius of the hot flow to a- m at the inner radius. Note that 
the rms(_E) dependence for the case of constant a would be 
constant as well. 

We then need to specify the fractional variability as a 
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0.01 



Figure 1. The function \Hq(x)\. The dotted and dashed curves 
give the approximations for low and high, respectively, values of 
x, given in Section ITl 



function of the radius. If we consider a perturbation with 
a given frequency, /, propagating with a constant speed, 
c s , towards a sink at the center (or from a central source 
outside), the wave amplitude will increase with decreasing 
radius as (see Vaughan & Nowak 1997; Nowak et al. 1999; 
Misra 2000), 



r(R) 



\H (2nR/\)\ 



(22) 



mK \H (2nB bl /X)\' 

where A = c s / f is the wavelength, Hq is the Hankel function 
of order zero, which modulus is given by 



\H (x)\ 2 = J (xf + Y ( X y 



(23) 



and Jo and Yo are the Bessel functions of order zero of 
the first and second kind, respectively. In writing equation 
12211 . we have assumed that the local rms level, r(R), is 
proportional to the wave amplitude, and normalized it to 
its maximum of r max at the inner disk radius. The func- 
tion | Ho (a;) | is plotted in Fig. Its argument is, e.g., 
unity for / = 3 Hz, c 8 = 0.01c, and R ~ 30GM/c 2 for 
M — IOMq. At x — > 0, it is logarithmically divergent, 
\H (x)\ ~ {1 + (2/tt) 2 [ 7 + ln(:r/2)] 2 } 1/2 , and for x > 1, 
\H {x)\ ~ (7tx/2)~ 1/2 , see Fig. Q Here 7 ~ 0.577216 is Eu- 
ler's constant. 

The above dependence of the wave amplitude as R^ 1 ^ 2 
can be understood as compression of the energy density 
of a constant-speed cylindrical wave with decreasing ra- 
dius, oc This increase becomes smaller when the wave- 
length becomes comparable with the radius, see Fig. We 
note that the sound speed in actual accretion discs gener- 
ally increases with decreasing R, which would flatten the 
dependence of equation I22H . A wave travelling through 
a geometrically-thin disc has c s oc R~ 9 ^ 20 and i? -3 ^ 16 in 
the gas-pressure and radiation-pressure dominated cases, re- 
spectively (e.g., Svensson & Zdziarski 1994). This would lead 
to r oc r~ 13 / 32 anc i _R -11 / 40 , respectively. The dependence 
for the gas-pressure dominated case is modified close to the 
inner radius for the case of the zero-torque boundary, lead- 
ing to c s — > at R —> Ri n , which would strongly amplify the 
wave. Furthermore, wave reflection from the inner boundary 
is neglected in the dependence of equation 1221 . A treatment 
of those effects is beyond the scope of this paper. 

The rms radial profile given by equation I22H has been 
obtained for a single frequency. Thus, that rms profile may 
apply for QPOs. However, the total rms dependence in some 
states follows that of the low-frequency QPO (Zdziarski et 
al. 2005), and equation 1221 is likely to be applicable more 
generally. We also note that perturbations propagating via 



waves will maintain coherence between emission at different 
radii (Vaughan & Nowak 1997), and thus this radial rms 
profile applies to the coherent variability only. 

On the other hand, the local fractional variability may 
be due to superposition of fluctuations (correlated or not) 
occuring at a number of radii (Lyubarskii 1997). In that 
case, it will be likely of a power-law form. We assume it 
here to depend as a power law of either the radius or the 
dissipation rate, 



r(R) = r max (R/R in )- p 
r(R) = r max [D(R)f /3 , 



(24) 
(25) 



where f3 is the index of the dependence. Note that equation 
125t with D(R) of equation 11611 is identical to equation 12 H . 
The profile of equation I22H at 2nR/X > 1 corresponds to 
f3 = 1/2 of equation 1241 1. Hereafter, we will consider only 
the values of f3 > 0. 



5 ANALYTICAL SOLUTIONS 

In the case of blackbody and Wien (see Appendix A) spec- 
tra with the torqued dissipation profile, equation I16L and 
Rovit — 00 we can obtain a number of analytical solutions 
for the coherent a, equation @. First, 



8tt 3 R 2 e - In (exp e - 

CTcoh = -r-Atj Tmax/iin 



(26) 



in the case of blackbody emission and /3 = 5/4. Here e = 
E/kT mm . For blackbody emission at a general value of /3, 
we find the limiting forms, 



Ccoh 



4/3/3 H 3 / \ 3 / 

C(8/3)r(8/3) 
e 5 / 3 ln(l/€) 
C(8/3)r(8/3) ' 
3e 5 / 3 

(4/3-5)<(8/3)r(8/3) ' 
1, 



£ < 1, < f 

e < 1, P = I 

e < 1, P > I 
e > 1, 



(27) 



where £ is the Riemann zeta function and V is the Gamma 
function. Note the power law index of 5/3 at low energies 
appears independently of the value of /3 as long at it is > 5/4, 
whereas it is 4/3/3 for (3 < 5/4. 

Then, the corresponding results for a Wien disc are, 

8-4/3 



a coh = ^-Ar^R 2 n E 3 e^-^ 3 



(8-40 \ 



and 

Ocoh 



4/3/3; 



' 8-4/3 



c 4/3/3 



V3 

r(8/3) 



e < 1; 
e»l. 



(28) 

(29) 
(30) 



Now, the low-energy power law is generally e 413 ^ 3 , in contrast 
to the blackbody case, where it appears only for < 5/4. 

In analogous way, we may calculate the dependences for 
the incoherent a using the integral approximations of either 
equation llOU or 11311 . For example, the Wien emission at a 
logarithmic coherence length yields, 

r l/2 f 16-8/3 „ \ 

i (AlnR) 1/2 (2e)^ /3 - \ 3 ' > . (31) 



(Tine 



2 ll/3' 
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The £< 1 limit is obtained trivially by replacing the incom- 
plete Gamma function by the complete one, which yields the 
same power law index as the corresponding coherent formula 
13011 , The expression is not valid for e > 1, where the ac- 
tual (Ji nc /F — > r max . This limit also applies to the blackbody 
emission, which case has then the e< 1 limit of 



~r max (AlnJ?)" 



1 4/3 i 

3 2 e -rr 2 | 



C 



13-8/3 x 



2r(8/3)C(8/3) 



0]' 



3e3 

"3 1 ! 

27 (4/3-5) Sr(8/3)C(8/3) 



/3< f, 



(32) 



It again has the same power law indices as the corresponding 
coherent cases, equation 1271 . Finally, /3 = 5/4 yields, 



47T -a o 1 

ff ta c = — ^£> max i?f„(Alni?)s 



ln(exp e — 1) 



e(exp e - 



- 1 



1/2 



(33) 



valid at any e except e ^> 1. At e <C 1, the factor in the 
second line approaches e _1 . 



6 VALIDITY OF APPROXIMATIONS 

Fig. |3 compares the accuracy of various approximations of 
Section |3J for blackbody emission with the local rms profile 
of equation 12 It with /3 = 3 and the temperature profile of 
equation 1161 . i.e., with the torqued inner boundary. Here- 
after, r max = 1, but the rms profiles below scale linearly with 
it. Also, in order to reduce the number of free parameters, we 
hereafter assume i? ou t = oo for blackbody and Wien emis- 
sion, which assumption has a negligible effect on the shown 
profiles as long as R out 2> 10 2 7?i n - 

The solid curve and the short dashes correspond to 
the case with full coherence, equations © and JSJl, and 
its discrete sum approximation by equations 1151 and 1141 
for Ami? = 0.1, respectively. We see that the approxi- 
mation begins to yield a co h/F lower than the integral for 
E/kTmax > 10. This is due to the discrete sampling of r, 
which results in the maximum sampled r being < 1 (in 
the middle of the innermost radial interval). The short/long 
dashes show the coherent rms in the case of local Wien emis- 
sion, equation 1291 . 

The dots correspond to the incoherent variability given 
by the sum of integrals, equation 1111 . at Ami? = 0.1, and 
the long dashes, to its approximation by the discrete sums, 
equations 1121 . 1141 . We see here the analogous effect to the 
one pointed above, with the discrete approximation sam- 
pling only r < 1, resulting in the rms lower than that due to 
the sum of integrals for i5/fcT max > 10. Then, the dots/long 
dashes show the rms approximation by equation 1131 with 
respect to the flux of equation JSJ. We see it starts to diverge 
for E/kT max > 10, which is due to using the integral to ap- 
proximate the discrete sum 1121 . Still, at E/kT max < 10 (for 
Aln 7? = 0.1), all the methods of calculating the rms yield al- 
most identical results in either the coherent or the inchorent 
case. Finally, the dots/short dashes correspond to the case 
with the linear coherence length with A(R/R[ n ) = 0.1 and 
equations which yields almost identical values of the 

rms as the corresponding logarithmic coherence length. 




100 



E/kT 



Figure 2. Comparison of various approximations to rms profiles 
for blackbody emission of a disc with the torqued inner boundary 
and a power law radial rms profile with /3 = 3. The solid curve, 
short/long dashes and short dashes show the case with full coher- 
ence, its Wien approximation, and approximation by a discrete 
sum with Alni? = 0.1, respectively. The dots correspond to the 
incoherent sum of integrals and Aln R = 0.1, and the long dashes, 
its approximation by the discrete sum. The dots/long dashes give 
in turn the approximation of the latter rms by an integral. The 
dots/short dashes correspond to the discrete sum with the linear 
coherence length with A(R/Ri n ) = 0.1. Here and in Figs. IHIS1 
the dissipation profile with the torqued inner boundary is used. 



7 EXAMPLE RMS ENERGY PROFILES 

Hereafter, we use the fully accurate formulae, i.e., with inte- 
grals and their sums. Equations @ and JJJ are used in the 
coherent case, and equations 1111 and JjjJ in the incoherent 
case. Subsequent figures are shown in the logarithmic scale 
of a/F, in which case the shape of a dependence is invariant 
with respect to changes of the normalization, r max . Also, it 
allows us to see the shapes of the dependencies at low values 
of a/F. 

For a given rms energy dependence, we have to choose 
the inner boundary condition (specifying the dissipation pro- 
file), equations 1161 or 1171 . and the radial rms profile, equa- 
tions 1221 (characterized by 27tRi n /A), or either 1241 or 1251 
(characterized by 0). For the power-law emissivity, equation 
1201 . we need to further specify Rout/Rin, c*in, « out , and 
E c . In the cases of incoherent variability, an additional free 
parameter is the coherence length, Aln J? (in the logarithmic 
case, which we will assume hereafter). 



7.1 Blackbody and Wien spectra 

We first consider the rms radial profile corresponding to disc 
waves, equation 1221 . For a given dissipation profile, the 
free parameter is the ratio of the inner disc radius to the 
wavelength, 2nRi n /X, dependence on which disappears for 
2nRi n /X > 1 as r(R) is then a power law with the constant 
index of —1/2 (Section^. Fig. [3] shows the coherent rms 
profiles for torqued and zero-torque inner boundaries and 
2nR in /\ = 0.01 and 1. 

Note that Figs. |5J|S| show the values of the rms up to 
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Figure 3. Example rms energy dependencies for the radial rms 
profile corresponding to disc waves (maintaining coherence). The 
short dashes and the solid curve correspond to the torqued in- 
ner boundary at 27ti?i n /A = 0.01 and 1, respectively, while the 
dots and long dashes show the corresponding results for the zero- 
torque boundary. The dependencies for that ratio > 1 are virtu- 
ally identical to those for = 1. 



E — lOOfcTmax, in order to show the shape of the depen- 
dence including their asymptotic behaviour. However, the 
blackbody emission becomes very weak for E > 15fcT max , 
and then other spectral components, e.g., Comptonization 
("Section 17.211 may dominate. To allow direct comparison of 
the shape, Figs. I^HHJhave the same range of a/F. 

Fig-Elshows a comparison between results in the black- 
body case with the power-law radial rms dependence and 
the torqued inner boundary for two values of f3 for both 
the coherent and incoherent variability. The rms(B) profiles 
change much more than in the case of the wave emissivity 
(Fig. |3J . The rms decreases with the increasing /3, and the 
incoherent variability is weaker than the coherent one for a 
given f3. Note that the rms at f3 > 5/4 has an asymptotic 
power- law form with the index of 5/3 at low energies, see 
Section IS] 

Fig. |S] compares the coherent rms in the case of local 
Wien emission, equation 12911 . with that for the local black- 
body for j3 = 2 and the torqued inner boundary. The black- 
body and Wien low-energy power-law indices are 5/3 and 
4/9/3 = 8/3, respectively, see see Section |5] 

Fig. |S] also shows the dependence of incoherent rms(_E) 
on Ami?. We see the rms increases with the increasing 
Ami?, and it approaches the coherent dependence at a a 
given -B/fcT max for a large enough Alni?. This is due to the 
incoherent emission being then dominated by the innermost 
disc element. 

Then, Fig.HJshows effects of the inner boundary condi- 
tion. We see that if we assume the zero-torque inner bound- 
ary but keep the rms(i?) as a power law, the resuling a/F 
are reduced by about an order of magnitude with respect 
to the case with the torqued inner boundary. On the other 
hand, assuming that the local rms does not directly follow 
the radius but instead the dissipation rate, equation 12511 . 
leads to the a/F approximately returning to the original 




E/kT m „ 

Figure 4. Dependence of the rms on j3. The solid curve and 
dots show the coherent dependence for /3 = 1 and 3, respectively. 
The short and long dashes show the corresponding incoherent 
dependencies for Alni? = 0.1. 




Figure 5. The coherent rms for the blackbody (solid curve) and 
the Wien spectrum (long dashes) for /3 = 2, and the dependence 
of the corresponding incoherent rms on Alni?. with the short 
dashes and dots showing the cases for Alni? = In 1.1 and In 1.5, 
respectively. 

values for torqued boundary condtion. Still, the incoherent 
rms is now lower by a factor of ~2 at high energies than that 
for the torqued boundary case due to the radial temperature 
profile being now less peaked. 

7.2 Thermal Comptonization 

We then consider the case of the e-folded power law, ap- 
proximating thermal Comptonization. Fig. Q shows the rms 
dependencies for f3 — 4 and various assumptions regarding 
the coherence and the profiles of dissipation and the local 
rms. The overall rms level strongly depends on those as- 
sumptions. All the dependencies are increasing, reflecting 
the assumption that the spectrum hardens with decreasing 
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Figure 6. Dependences of rms profiles on the disc inner boundary 
condition. The solid curve shows the coherent dependence for 
fi = 2 for the dissipation profile with the torqued inner boundary. 
Then the dots show the corresponding profile for the zero-torque 
inner boundary and the local rms being a power law in radius, 
equation 1241 . The dots/short dashes show the analogous profile, 
but with the local rms following the dissipation rate, equation 
1251 . The short dashes, long dashes and dots/long dashes show 
the corresponding incoherent profiles for Alni? = 0.1. 

radius, though much flatter than in the blackbody case. We 
can obtain the rms increasing with radius for either the spec- 
trum softening or the local variability decreasing with the 
decreasing radius. 

Fig. |S| show the effect of varying f3 for both coherent 
variability and incoherent one with Ami? = In 2. The in- 
creasing (5 results in decreasing rms. 



8 APPLICATION TO GRS 1915+105 

We compare here our theoretical models with local black- 
body emission to the rms from an observation of the micro- 
quasar GRS 1915+105 performed by the RXTE/PCA on 
1999 April 21. Its variability properties have been studied 
in detail in the Zdziarski et al. (2005). It belongs to the ui 
(Klein- Wolt et al. 2002) variability class. It is characterized 
by a very bright, ultrasoft spectrum, which is well fitted by a 
disc blackbody Comptonized by a relatively cold, Thomson- 
thick plasma (Zdziarski et al. 2001). Still, we have found 
that the 3-20 keV spectrum can be approximated within 10 
percent by an absorbed disc blackbody alone with the max- 
imum temperature of fcT max ~ 2.3 keV. Thus, we apply to 
it the formalism developed above for such spectra. Details 
of determination of the rms profile are given in Zdziarski et 
al. (2005). Fig.Elshows the 3-18 keV rms, obtained directly 
from the lightcurve with 1/128 s time bins. 

We have performed an extensive search of theoretical 
models fitting the data. We have found, in particular, that 
we can rule out models with the rms profile corresponding 
to the amplitude of a disc wave, equation 1221 . The reason 
for it is its rms(_E) dependence being much too flat (Fig. 

than the observed dependence. The data require, in fact, 
an rms profile much steeper than j3 = 1/2 of the wave pro- 



T 1 1 1 — I I I I | 1 1 1 1 I I I I | 1 1 1 1 I I I I 
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0.1 1 10 100 
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Figure 7. Example rms dependencies for an e-folded power law 
with the spectral index changing from a ut = 0.9 at 100-Ri n to 
o in = 0.4 at _R in for E = 0.1 keV, E c = 100 keV, and /3 = 4. The 
solid curve, dots and short dashes correspond to the power-law ra- 
dial rms profile and the torqued inner boundary, the same profile 
for the zero-torque boundary, and the the rms profile following the 
latter dissipation profile, respectively. The long dashes, dots/short 
dashes and dots/long dashes show the corresponding dependen- 
cies, respectively, for incoherent variability with Ami? = 0.1. 




E [keV] 

Figure 8. Example rms profiles for an e-folded power law with 
c»in, a ut, Rout/Rin, Eq, E c as in Fig. Q the dissipation pro- 
file with the zero-torque inner boundary and the rms governed 
by its dissipation profile. The solid curve and dots correspond 
to coherent variability at /3 = 1 and 3, respectively. The short 
and long dashes show the corresponding incoherent variability at 
AlnR = In 2. 

file. Our best model, shown in Fig- El h as P = 4, incoherent 
variability with Alni? = 0.1 and the dissipation and rms 
profiles of equations 1171 and 1251 . respectively. The flat- 
tening seen at low energies then requires an addition (in 
quadrature) of a variability pattern with the rms indepen- 
dent of energy (corresponding to a varying normalization of 
the entire spectrum) at the amplitude of cr cons t/F = 0.115. 
We stress we have approximated here the photon spec- 
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Figure 9. The rms observed from GRS 1915+105 in an ultrasoft, 
blackbody-dominated, state. The solid curve show the theoretical 
model with incoherent variability with Alnii = 0.1 in a flow with 
a zero-torque inner boundary, the radial power-law rms profile 
following the dissipation rate at j3 = 4 and normalized to r max = 
1.35. A weak rms component constant with energy has also been 
added (see text). 

trum as a disc blackbody, which only roughly approximates 
the actual spectrum. Taking into account this inaccuracy 
will modify somewhat the model rms dependencies. There- 
fore, we we have not systematically fitted the rms data, e.g., 
by the least-square method, as the resulting small fit im- 
provements may be less significant than the inaccuracies in- 
troduced by our assumption above. 



9 DISCUSSION AND CONCLUSIONS 

We have presented a formalism to calculate rms spectra from 
accretion flows under the assumptions that each radius is 
characterized by a given value of the local rms and a given 
spectral form (Section |3J. We do not constrain a priori the 
profiles of the local rms, but it is likely that the flow be- 
comes more turbulent and unstable with decreasing radius, 
which results in the rms increasing with decreasing R (at 
least down to some radius). In Section 0] we consider three 
specific radial emissivity profiles. In the case of a local black- 
body or Wien spectra, the peak energy increases (at least 
for most radii) with the decreasing radius, which leads to 
rms spectra increasing with energy, see Section \7. II This is 
indeed observed in an example of a blackbody-like emission 
we consider (Section|HJl. For the case of local Wien emission, 
we obtain analytical formulae for the integrated disc spec- 
trum, which we find to be a good approximation to the disc 
blackbody spectrum (Appendix A), and the rms, Section [H] 
For both blackbody and Wien spectra, we obtain formulae 
connecting the power-law index, (3, of the radial rms profile 
with the low-energy index of the energy rms dependence, 
equal to either 4/3/3 or 5/3 (Section^. Thus, measurements 
of the indices of rms energy dependence can constrain the 
variability properties of the emitting accretion disc. 

We have also considered the local emission due to ther- 



mal Comptonization (Section |1J . If the seed photons are 
emitted by an outside cold disk (as is likely to be the case 
of the hard state of black-hole binaries), the power-law part 
of the local Comptonization spectrum will harden with the 
decreasing radius. Such a behaviour is indeed observed in 
some sources (Revnivtsev et al. 1999). Again, integration 
over radii yields rms spectra increasing with energy (Sec- 
tion 17.21 , though with an overall dependence much flatter 
than in the blackbody case. Such a behaviour is also ob- 
tained in models considering in more detail propagation of 
disturbances towards the central black hole (Zycki 2003; P. 
T. Zycki, personal communication). 

On the other hand, some black-hole binaries in the hard 
state show X-ray rms spectra decreasing with energy (e.g., 
Wardzihski et al. 2002). This is likely to be caused by varying 
shape of the spectrum emitted at a given radius (Section 0, 
which effect is not included in the formalism presented here. 

We have also not considered variability related to the 
orbital motion at a given radius. This effect may be impor- 
tant observationally in AGNs (Turner et al. 2004), where the 
orbital motion can be well resolved by X-ray observations. 

Given the theoretical uncertainty about the level of tur- 
bulence in accretion discs as a function of radius and the 
degree of correlation between emission from different radii, 
our models rely on a number of assumptions and free pa- 
rameters, e.g., about the radial profile of the rms, see Sec- 
tion |7| Still, data of sufficient quality can strongly constrain 
the models. For example, blackbody rms energy dependence 
flattens at a much lower energy in the coherent case than in 
the incoherent one (Section 17. 111 . 

We stress that rms studies can, and should, be com- 
bined with other observational constraints. The components 
present in time-averaged energy spectra can be determined 
by spectral fits. Furthermore, the degree of coherence be- 
tween emission at different energy bands (and, presumably, 
different radii), can be found by cross-correlation of the cor- 
responding lightcurves. A more detailed diagnostic is the 
coherence function (Vaughan & Nowak 1997), the cross- 
correlation coefficient as a function of Fourier frequency. 
Also, coherence between spectral components can be tested 
directly from the rms energy dependence, as the rms(i5) 
from coherently varying components would be simply the 
weighted average of the rms of the components (Section 
On the other hand, two components varying independently 
would give a dip in the rms(i?) at the energy where their 
spectra intersect, due to the partial cancellation of their in- 
dependent variability. Such dips appear, e.g., to be absent 
in the x states of GRS 1915+105 (Zdziarski et al. 2005). 

In future work, the various approaches to calculating 
observed rms spectra can be combined. E.g., more than one 
component with radial dependencies of rms (Section |5J can 
be considered (Zdziarski et al. 2005), or the radial depen- 
dence can be combined with time evolution at a given radius. 
Then, the rms from time evolution of a uniform Comptoniza- 
tion source due to variability of the plasma parameters is 
considered by Gierlihski & Zdziarski (in preparation). 
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Figure Al. The disc blackbody spectrum (solid curve) and the 
disc Wien spectrum (dashes) for the torqued inner boundary (in 
the EFe representation). 
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APPENDIX A: DISC BLACKBODY AND DISC 
WIEN SPECTRA 

In the case of the torqued inner boundary, equation 11611 . 
and .Rout = oo, the disc blackbody spectrum (equations 
ED) is, 

Jl Bxp(et)-H 

e-»/ 3 r(|)c(|), e<l; ( A1 ) 

e _1 exp(-e), e > 1, 

where e = E/kT mBix , £(8/3) ~ 1.28419, and L(8/3) ~ 
1.50458. 

On the other hand, the blackbody spectrum in the Wien 
limit becomes, 

F(R) = AE 3 exp[-E/kT(R)]. (A2) 

This is actually the equilibrium photon distribution reached 
in a purely scattering medium (without stimulated scatter- 
ing) with infinite optical depth, and thus it may be a bet- 
ter approximation to emission of effectively-thin (but opti- 
cally thick to scattering) accretion discs (Shakura & Sunyaev 
1973) than the blackbody distribution. Wien-emitting discs 
have been considered, e.g., by Misra, Chitnis & Melia (1998) 
and Misra (2000). 

With the same assumptions as for equation HAIL the 
disc Wien spectrum can be expressed in a closed form, 

F=^^V 8 / 3 r(|, e ), (A3) 

Notably, this spectrum is very close to the blackbody spec- 
trum, see Fig. IA1I Its limiting forms are the same as those 
in equation lAH except for no ("(8/3) factor at e <C 1. 
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